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Abstract 

Efimov physics relates to 3-body systems with large 2-body scattering lengths a s and small effective 
ranges r s . For many systems in nature the assumption of a small effective range is not valid. The present 
report shows binding energies E3 of three identical bosons calculated with 2-body potentials that are fitted 
to scattering data and momentum cut-offs A by inverse scattering. Results agree with previous works in the 
case of r s < o s . 

While energies diverge with A for r s = 0, they converge for r s > when A >~ 10/r s . With a" 1 = the 
converged energies are given by E% — C^r~ 2 with n labelling the energy-branch and calculated values 
C^ 0) = 0.77, C^ 1 ' = .0028. This gives a ratio ~ 278 thus differing from the value ~ 515 in the Efimov case. 

Efimov's angular dependent function is calculated. Good agreement with previous works is obtained for 
r s <C a s . With the increased values of r s the shallow states still appear Efimov-like. For deeper states the 
angular dependence differs but is independent of r s . 

1 Introduction 

Systems involving particles with 2-body interactions at or close to the Unitary limit have become of specific 
interest in the physics community during the last several years for reasons that have repeatedly been pointed out 
in numerous publications related to both atomic and sub-atomic problems. The Unitary limit is here defined 
as being that for which the scattering length a s and the effective range r s is infinite and zero respectively. 

In this regard the bound system of three bosons is of special interest. This was first brought to the attention 
by the works of Thomas [5] and of Vitali Efimov [3J 0]. Bosons interacting with a resonance in the 2-body state 
(i.e. aj 1 ~ 0) result in a strongly bound three-body system and with a spectrum of loosely bound excited 
states. 

The literature has brought to the attention several systems in nature for which these findings are relevant. 

On the theoretical side the application of Effective Field Theory (EFT) methods have been proven very 
powerful^ H3 |H1 HI [ID] for a theoretical interpretation. The main focus has been on short-ranged potentials 
at or near the Unitary limit for which r s <C a s , but more recent EFT publications include effective range 
corrections. [H] [12] 

A rather different approach is used here. The inverse scattering formalism for separable potentials is used to 
construct potentials as functions of the scattering parameters as well as 2-body binding energy, and renormalised 
by a cut-off in momentum-space. The potential in the Unitary limit, a renormalised ^-function (in coordinate) 
space, is of specific interest. The results presented below agree with previous Efimov physics calculations for 
r s = 0.. The r s > results do in general show some qualitative differences. 

The present investigation focuses on Efimov energies and the Efimov function, dependence on scattering 
length and in particular on the effective range. Some preliminary results were given in ref.|13j 

Scaling properties are a main thread in this work. 

In Section 2 is found a presentation of the necessary tools which are the Faddeev equation and the inverse 
scattering method. Section 3 show results of numerical calculations in 4 subsections with 8 Figures. Section 4 
is a summary and some discussion of the results. Some relations regarding off-shell scattering and three-body 
forces are shown in Appendix A, while Appendix B shows an important relation pertaining to scattering and 
the separable interaction in the Unitary limit. 



2 Formalism 

2.1 Two-body Separable Interaction 

The input for the calculations in this report are scattering phase-shifts S(k) related to the scattering length a s 
and an effective range r s by 

fccot<5(fc) = -— + \r s k 2 (1) 

Qj 8 Z 
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A momentum cut-off A will also be defined so that 



*(* > A) = 

With the scattering parameters chosen so that the phase-shifts do not change sign, as will be the case here, 
they can be reproduced by a rank-1 separable potential 

V(k,k') = -v(k)v(k') (2) 
that is obtained by inverse scattering from[l4] 

4-7T 

v 2 (k) = — sin 5(k)\D(k 2 )\ (3) 



where 



^ 2 ) = ^^exp 



2 f A 

n ' L k 2 -k' 2 



(4) 



where V denotes the principal value. The interaction is fully defined by the phase-shifts and the two-body 
binding energy Eb- Note that the two fits, to scattering and to binding energy are independent. The equation 



sfE^ = (1 - yl - 2r s aJ 1 )- (5) 

relates Eb to [a s ,r s ] but only for r s <C a s and it is not useful for the calculations below. In the limit r s = it 
reduces to 

E B = a; 2 

which will be used also when r s ^ together with Eb — for a s < 0. This simplifies the comparison with the 
r s — calculations in these preliminary calculations. The major importance is that the 2-body system has a 
bound state. 

When dealing with some specific physical system one (in general) knows or assumes some binding-energy, 
but this is not the case here. The expression ([5]) for the binding energy does not introduce another parameter 
in the theory, being a function of a s and r s only, but if in some specific case one would tune the binding enrgy 
say by some shape-parameter there would be another parameter to consider. 

2.1.1 Unitary Limit 

A well-known reason for assuming a separable interaction in 3- and other many-body works is that it in general 
simplifies the formalism as well as the numerical calculations. There is however a particular case where the 
interaction IS separable. Namely that when the 2-body system has a bound state at or close to zero [IB], which 
is the case in the Unitary limit. 

The principal value integration in eq. @ can be done analytically for <!>(£;)= const ant. The Unitary limit is 
a special case with 5(k) = ~, for which the rank-1 separable potential v u (k) is given by [TS] 

vl(k) = - f ** (6) 

(A 2 - fc 2 )2 

If A 3> k one finds 

vl(k) -+ ~ (7) 

In this limit, but only in this limit , the unitary interaction is independent of momentum and a 5- function in 
coordinate space with the strength inversely proportional to the cut-off. But for the finite values of A needed 
for computations, renormalisation is required resulting in v u as given by eq. ([B]). Note that there is then an 
abrupt increase in strength and a singularity as k — ¥ A. This is required to preserve the condition 8 = ? for all 
k < A, as required for the Unitary limit. (See Appendix B.) 

The singularity can cause numerical problems that however can be largely overcome by proper computational 
methods. The substitution 

k = Asin^, 



2 Calculations showed no qualitative difference if using eq. f5l . 
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Figure 1: Lowest curve (red on-line) shows the diagonal of the Unitary interaction in momentum-space. The 
two other curves are, from below, the potentials for a s = 10 (green on-line) and for a s = 20 (blue on-line). All 
three are with the effective range r s — 0. 



with 9 the new variable, can for example be helpful yielding 

which is used successfully in the computations. 0. Fig. [T] shows the Unitary interaction (lowest curve) and for 
comparison two other potentials with a s — 10 and 20 respectively. 



2.2 The Many-body System 

A rank-1 potential is often sufficient to fit the on-shell data at low energy as is done in this report. The off-shell 
and off-diagonal parts of the T- matrix is then fixed by default. 

On-shell relates to the asymptotic form of the two-body scattering wave-function, while off-shell relates to 
the interior wave-function, i.e. to correlations in the two-body system. One may however like to change the 
off-shell behaviour. This can be done by extending to a rank-2 (or higher) potential (which implies introducing 
another parameter), while preserving the on-shell fit. This provides a practical tool for exploring off-shell effects 
on many-body properties. □ 

The 3-boson system collapses in the Unitary limit (the energy diverges) with A. This is well-known. In 
ref.[5] this collapse was remedied by adding a 3-body force. It was shown in ref.[13] that the same could be 
accomplished by a change in off-shell scattering with a rank-2 potential as described above. 

In each of these methods the high energy (ultraviolet) region of the interaction was modified; in one case by 
a three-body force in the other by a change in off-shell scattering. 

Although the two methods have a similar 'end-result' it is important to understand that in each case there 
is a different physics involved. (Appendix A.) 

Findings below show that the effective range also has the effect of preventing a collapse. This is not 
completely un-expected as it explicitly introduces a finite range parameter into the system. 

3 This is in particular useful because dk — > A cos 8d8 which can eliminate the cos0 in the denominator of JHJ 
4 It was for example used by the author in earlier work on the triton. |17| 
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2.3 Faddeev Equation 

The Faddeev three-boson equation for a spin- independent rank-1 separable potential is given by 



2 /■ 2A «(|k + lq|H|q+lk |) 

l(E 3 ~lqi)J q 2 + < l -k + k 2 -E 3 



A 



with 

^(«) = 1 + 77^7 / ^(fclfs-fc 2 )- 1 /^ 2 ^ (10) 

2tt 2 Jo 

With the phase-shifts defined, as announced above, in terms of two scattering parameters and the cut-off A 
one will have E 3 — E 3 (a s ,r s , A). In the present investigation the scattering parameters [a s ,r s ] are considered 
internal parameters defining the 3-body system, while [A] is considered an external parameter if possible chosen 
larger than the maximum momenta of the bound particles. This means in general that one has to choose 
A > 1 /R with R being the physical size of the system. This size is a 'functional' of the interaction (the internal 
parameters) and so is therefore A. As an example, the three-body system collapses in the Unitary limit, with 
A — > oo , a well-known divergence. 

Three separate cases will be considered here with results presented in three separate sections. The first will 
be the Unitary limit i.e. both r s and a^ 1 =0. The second is aj 1 ^ while r s = and the third the more 
general case, both r s and aj 1 ^ 0. The purpose of the present investigation is in particular to investigate the 
latter case, with r s > 0. 



3 Numerical Results 
3.1 Unitary Limit 

The Unitary limit is a special case. After a change of variables (see above), k — > A sin 6, the function X(s), eq. 
(ITU1) . becomes 

X(s) = l + -V { 2 sin 2 6{s- sin 2 O^dO (11) 
t Jo 

The integral, done analytically for s < yields 

X(s) = —^= (12) 
Vs — s 

(See Appendix B for s > 0). The only free parameter in the Unitary limiit is A so that one has E3 = Ea(A). It is 
then convenient in this case to choose momenta and energies in units of A and A 2 respectively by substitutions 
k -> fcA q -» qA and E 3 = E U A 2 in eq. (J9j> to get with eq. jH} and s = E u - f q 2 



2yf^~~s /- 2 ?; u (|k+iq|K(|q+ik|; 



^ - -^1— J 7 -i x{k)dk (13) 

(n) 

With n labelling a specific state with energy Eu the eigenvalue spectrum in this Unitary limit was calculated 
with the result (in units of A 2 ): 

E u °) = -.1325, = -.0002520, E u 2) = -.0000004568 

One finds 

E u ^/E u ^ = 515 

agreeing with the Efimov result while 

E^/eW = 630 

This latter discrepancy is ascribed to computational inaccuracies that increases for the smaller energy. (E u = 
-.0000004893 would be the 'correct' value.) 
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Figure 2: Energy as a function of scattering length a s while the effective range is r s — 0. Following convention, 
K = k = v^3- The lowest curve (red on-line) is for n = 0. The two other curves (green and blue respectively 
on-line) are for n = 1 and 2 respectively. The crosses on the dimer-energy line, shown by the broken line, are 
explained in the text. A = 10 in this plot. 

3.2 Non-zero Scattering Length but Zero Effective Range 

In this second scenario with r s — while a s ^ one will have E 3 — E s (a s ,A). This case has been the focus of 
numerous investigations since the early works of Thomas and Efimov; in particular as relates to experimental 

studies of cold atoms and the ability to tune the scattering length with the aid of Feshbach resonances. 

i _ i 

Results of Faddeev calculations are shown in Fig. [2] Following convention, the figure shows E 3 4 vs a s 2 
(with appropriate signs) for three separate branches cutting the a s = axis at the values shown above. With 

k* = Es(a^ 1 = 0) for each branch it has been found, e.g. [7 , that the energys should approach zero for 
a' g ~ —1.56k* 1 . These scattering lengths are indicated by a "Y" for each of the states and largely agree with the 
calculations. Furthermore, the trimer and dimer energies are estimated to coincide at a* = 0.0707645086901k~ 1 - 
These points are shown by crosses in Fig. [2] for the two shallow states. For the lowest state the corresponding 
point is predicted to be at a J 1 as 51. 

It is convenient to choose A as the unit of momentum and introducing the dimensionless variable a s A. The 
phase-shifts and therefore the potential are then functions of this variable. So in addition to the A 2 factor the 
energy will also be a function of the same variable. One will have 

E { 3 n) = A 2 F^(1/Aa s ) (14) 

where as before n labels a specific (Efimov) branch. (The reason for choosing the inverse of Aa s becomes obvious 
below.) While Fig. [2] shows the energy vs a s for A = 10 the energy for any other value of A and a s can be 
obtained from eq. ([M)) . Of interest is then the function F^- n \ The full line in Fig. |3] shows sj F( n )(l/Aa s ) 
for n = 0, the deeper state in Fig. [21 It is seen that this function is to some approximation a linear function 
with the slope (derivative) being approximately one. The function F^ is a (nearly) scaled copy of the 
scaling factor being 22.7, as shown by the dotted line in Fig. [31 The reason for displaying these functions is for 
comparison with the similar situation when r s ^ that is shown in Sect. 3.4. 

From the functions one can calculate the energies E$ for any other variables a s and A for each branch 
n. These functions can of course also be calculated directly from the Efimov Universal function, A(£). It is 
related to the energy E 3 by 
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Figure 3: The function V fW is shown for the two deepest states. The full curve (red on-line) is for n — while 
the broken (green on-line) is for n — 1. This curve has been scaled up by the Efimov factor 22.7 as a result of 
which the two curves are close. Note also the approximately linear function of l/a s A. 



E { 3 n) = -a-2 + Kl(e~^) n - no e^r (15) 

where sq w 1.00624 and n labels a specific Efimov branch and no that branch for which k 2 = E 3 (a^ 1 = 0) . 
The angle £ is defined by 

tan£ = —a s VE s ^ (16) 

It is esily seen that because of the the scaling relation shown above A(£) is indeed a function of £ only, 
independent of A. 

A(£) was calculated with twelve significant figures by Mohr[9]. Earlier results by Braaten et al [6] were 
shown parametrized in their publication. 

The dots and crosses in Fig. @]show A(£) calculated from eq. (|15p using the same energys as used for Fig. 
[2J The crosses refer to the energy E^ , the deepest state, and the dots to E^\ the next shallower state. The 
broken line is from the parametrisation given in ref [6] . The present calculations are not of high-precision but 
the difference between our result and those of ref. [B] are still less than a few percent. 

3.3 Non-Zero Effective Range and Zero Scattering Length 

After having verified that the r s — calculations with the separable interactions agree with previous works the 
rsesults with r s ^ can now be presented. 

The scenario changes drastically with the introduction of a non-zero effective range. As a preliminary to 
the more general case the scattering length is first assumed to be infinite. The remaining variables are then r s 
and A. While in the previously considered case E 3 oc A 2 associated with a collapse of the system (in coordinate 
space), the range-parameter r s now provides some finite size R of the system under consideration. As is to be 
expected computations of E3 as a function of A now show a convergence for A > 1/R. With A providing a 
momentum scale and following arguments above one can expect a relation of the following form 

^=A 2 F (n) (r.A) (17) 
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Figure 4: The function A(£) is shown. The crosses (red online) are for n 
n = 1. The broken line (black on-line) shows the results of ref.[6]. 



the dots (green on-line) are for 



In order for convergence one must then have 



^ ( " )(rsA) ^(W (18) 
for large A. The energy of the three body system will then scale with r s as 

E = 

' s 

for A large. The coefficients CW were calculated to get = 0.770 and = 0.00298. The ratio C^/C^ ~ 
278 differs from the related ratio of 515 for Efimov states. 

Fig. H] shows and 278 * E { 3 1} as a function of A for r s = 0.1. The energy for any other values of the 
parameters r s ^ and A can be calculated by the help of eq. (p~7|) . Calculations have shown that the energy 
converges to the asymptotic value at a critical value 



1 



A c ~ 10/r, 

as in FigJS] where one finds A c ~ 100 for r s = 0.1. 0. 



3.4 General Case 

The main inspiration for the present work has been to find the effect of a non-zero effective range r s > together 
with an inverse scattering length aj 1 ^ as opposed to the case when r s = 0. 

There is an additional factor to address here, the dimer (two-body) binding energy. In the inverse scattering 
theory with separable interactions it is chosen separate, independent of the fit to scattering data. Following the 
discussions in Sect. 2.1 the choice here is Eb = aj 2 - Any other choice might change our results but (probably) 
only in a qualitative way. In cases of real physical systems one would of course choose the appropriate dimer 
binding energy assumed known for the particular system under consideration. One consequence of this choice 

5 Fig. 3 in ref . [13] shows E vs A for r s = .0, .03, .05 and .1 with convergence at A c ~ 10/r s . The convergence A c , as a function 
of r s is also seen. The curves for a^ 1 = (excluding the r s = 0) are related by the scaling in eq. Ill7jl 
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Figure 5: The energies E^ (red on-line) and E^ (green on-line) are shown here as a function of the momentum 
cut-off A with the effective range being r s — 0.1 while aj 1 = 0. The energy E^' is multiplied with the factor 
"278" (see text) and the two energies practically overlap. The energy for any other [A, r s ] is obtained by the 
scaling shown in the text. Of particular interest is the asymptotic value, i.e. A > A c , that scales as r~ 2 . 



is that no additional parameter is introduced in the calculations. With A > A c (see above) the energy is 
then a function only of a s and r s . Choosing momenta in units of r" 1 the potential will be a function of the 
dimensionless parameter a s /r s . The energy £3 then has to have the form: 

BfUr-^^a- 1 ) (19) 

where the functions are to be determined computationally!^ The scaling shown here is analogous to that 
shown in eq. (|14l) . with r s replacing A and providing the new momentum scale. 

Fig. [3J showed functions F^(l/Aa s ). The similar functions F[ n \r s / a s ) are shown in Fig. [6l The function 
Ff 1 ' maps onto F' 1 ' in Fig. O with a scaling factor of ~ 3. So Efimov physics seems to apply for the shallow 
state. With regard to the deeper state there are some qualitative similarities but there are definitely quantitative 

differences. Rather than having a slope ~ 1 as in the Efimov case, Fig. |6]shows \l F^ to have a slope of ~ f .5. 

All the results in this section, for r s 7^ can be obtained by scaling from the functions F^ shown in Fig. 
[SJ As a first example is shown Fig. [JJ the energy vs a s for r s = 0.1. These curves are seen to be qualitatively 
similar to the 'Efimov-curves' shown by Fig. [2] From these data a function A'(£) was calculated using eq. 
(|T5J) with the result shown in Fig. [5] There is a factor (e~^) n ~"° = 515™""° in eq. (JTSJ) for A(£). For the 
purpose of calculating A'(£) it was replaced by 278™ _ ™°. The function F^ generated the result shown by the 
crosses, while F^ the dots. The dots, referring to the shallowest state spans r s /a s < 0.8. They coincide with 
the Efimov A(£) (the broken line), while there is a definite difference for the deeper state derived from F^°\ 

The situation here is however anologous to that in sect. 3.2. There, the consequence of the scaling expressed 
by eq. (|14l) had as a result that the calculation of A(£) was independent of A. The scaling expressed by eq. 
(|19p can formally be obtained by replacing A in eq. (|14l) by r s . The function A'(£) is similarly independent of 
r s - 

6 Like in Sect. 3.2 for the function F the argument of Fi is chosen as the inverse of a s /r s . 
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Figure 6: The full curve (red on-line) shows y F± as a function of r s a s 1 . The broken curve (green on-line) 

shows \J scaled up as in Fig. [3] It is seen that it (the broken) only spans values of r/a < 1 (note the 
scaling factor "22.7") so that Efimov physics would or might apply here. The full line, for n — 0, (red on-line) 
includes energies for r s /a s = 20 and there is as a consequence an apparent deviation from Efimov physics here. 
This is substantiated by the results shown in Fig. [H 




-7-1, , , , , , , , , , , J 

-b -2 -1 G 1 2 b 4 b 6 7 8 

V1 /a 

Figure 7: Energy (K) as a function of scattering length for r s = 0.1. The energy ratio between the two states 
at 1/a — is here ~ 278. This factor is independent of r s . It is approximate and improved calculations are 
called for. 



9 




Figure 8: The function A'(£) calculated from eq. (p~5|) . with "515" replaced by "278" as described in text, 
for each of the two states shown in Fig. [7J This function is independent of the value of the effective range r s . 
The dots (green on-line) refers to the shallow state, while the crosses (red on-line) the deeper. The broken line 
(black on-line) shows the Efimov function using the parametrization in ref.[6]. 



4 Summary 

Bound state energies E3 of the 3-boson system with 2-body interactions at or close to the Unitary limit were 
calculated with the Faddeev equation in momentum-space. Separable potentials were obtained by inverse 
scattering as a function of scattering length a s , effective range r s and momentum cut-off A. 

In the Unitary limit (r s = aj; 1 = 0) three branches (labelled n = 1 — 3) were identified with energies (in 
units of A 2 ), -.1325, -.000252, -.00000045668 respectively. The ratio E { 3 n ' 1} / E ( 3 n) ~ 515 in agreement with 
theory. 

With r s = 0,a s 7^ (Sect. 3.2) the energy for each branch scales as E 3 n ^ = A 2 E^ n \a s A). Efimov's function 
A(£) was calculated separately for two branches (n — 1,2) with the result shown in Fig. @] agreeing with 
previous works H] . 

While the energy oc A 2 for r s = the situation changes when r s > 0. One still finds the energy oc A 2 for 
small A but it converges for A > A c ~ 10/r s as seen in Fig. [5] (Sect. 3. 3) for a^ 1 = 0. It was shown that the 
energy converges to E^{n) — C*-™Vj 2 with ri ^ being the only remaining momentum-scale with A > A c . The 
ratio C^™ -1 * 1 /C^ is calculated to be ~ 278 differing from the corresponding Efimov value ~ 515. 

It was shown in Sect. 3.4 that if both r s > and aj 1 ^ the energy can be written as E 3 — r~ 2 F± (rgaj 1 ). 
The function was shown and a function A'(£) was calculated analogous to the Efimov function A(£). Fig. 
[5] shows it to overlap with the Efimov for the shallow state. This might be explained by the fact that for this 
state, r s is still < aj; 1 . For the deeper state, there is however a definite difference even for that range of a s . 
An important result following from the expression ([TT))) for the energy J5g (r s ,a s ) is however that A'(£) is 
independent of r s > also for these states. Further investigation is asked for here. 

There would (probably) be some quantative changes in the results above if a different algorithm were used 
for the binding energy as a function of a s and r s , but this all relates to a specific system, beyond the scope of 
the present investigation. It would in general introduce another parameter into the theory. 

The Efimov physics, valid for r s /a s <C 1 is well understood and explained by various methods. Particularly 
powerful is the EFT approach with recent and on-going research [lll 112] extending this method to finite values 
of the effective range. 
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The main purpose of this work is to provide an alternative approach to the the problem at hand and in 
particular to explore the domain of r s > 0. The results were obtained with rank-1 separable potentials. It is to 
be expected that any other potential fitted to the same scattering data would give the same, or very similar, 
results. 

Although much care was taken, the computing was not of 'high precision' and could be improved upon, but 
this is not expected to have an impact on qualitative conclusions. 



5 Appendix A 

The 1Z (or T)-matrix relates the potential to observable scattering data while the potential in itself is not 
an observable. The phase-shifts 8(h) are related to the reactance matrix 1Z and indirectly to the two-body 
interaction V by 

H{u) = V + VP^^KK (20) 
< k\K(w = k 2 )\k >= i tan(5(fc) 

As a consequence only diagonal elements of the 72.-matrix can be obtained from the phase-shifts and only for 
lu = k 2 . In addition, experimental information puts an upper limit on k. A large part of the matrix is therefore 
beyond direct experimental reach. Many-body theories do in general need more, requiring some other input, 
e.g. by potential models. 

An important example is off-shell information. Off-shell implies lu ^ k 2 . It is therefore relevant to consider 

^ = (vf-Sl* = (21) 

* = $ + l -rF,V^ 

lu — k' z 

where $ is the unperturbed wave-function. 

These relations were already used by Moszkowski and Scott. [TBI IE] I w is often referred to in the literature 
as the wound-integral. The 'correlated' wave-function ^ relates to details of the interaction (even beyond what 
may be available from experiments.) 

Assume a many body system with 2-body interactions only; no 3-body (or rt-body, n > 2) forces. Then 
consider two particles (i,j) interacting in this system. The interaction would be off-shell, lu ^ k 2 , mainly 
because of binding-effects, meaning that other particles (k) affect the the interaction between (i,j). So to 
calculate something, e.g the total energy, one would have to sum over (i,j,k) that 'looks' as if one is dealing 
with a 3-body force , although 2-body only was assumed. It is then more appropriate to say one is dealing 
with 3-body terms not forces . Whether the 3-body term is 'important' or not is another question. There 
are of course in principle also 4,5 etc body terms. That depends on the system in question including density, 
temperature etc. 

The above is to illustrate the popular statement of equivalence between the off-shell effect and that of a 
3-body force, even though they have a different origin. The first is a medium-effect, the second relates to internal 
degress of freedom of the particles. Either or both may be relevant for a specific system in nature. In either 
case there will for example be the same energy vs density functional. 



6 Appendix B 

The scattering phaseshift is related to the reactance matrix in Appendix A. For a rank-1 separable potential 
one finds 

< k\K{u = s)\k >= ^ (22) 
l[s) 

In the Unitary limit I(s) is given by eq. (| 12[) with s < in the Faddeev equation, eq. (|13[) . But for free space 
2-body scattering one has s = k 2 so that with momenta k < 1 (in units of A) one will have < s < 1. The 
integral, can again be done analytically and yields 

1(0 < s < 1) = 
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so that the Reactance matrix element 

< k\R\k >= - tan<S(fc) -> oo 
k 

i.e. S(k) = ^ for ALL momenta < k < 1 and the condition for a unitary interaction with a cut-off A in 
momentum-space is satisfied. 
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